Theory of Feshbach molecule formation in a dilute gas during a magnetic field ramp 
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Starting with coupled atom-molecule Boltzmann equations, we develop a simplihed model to 
understand molecule formation observed in recent experiments. Our theory predicts several key 
features: (f) the effective adiabatic rate constant is proportional to density; (2) in an adiabatic 
ramp, the dependence of molecular fraction on magnetic field resembles an error function whose 
width and centroid are related to the temperature; (3) the molecular production efficiency is a 
universal function of the initial phase space density, the specific form of which we derive for a 
classical gas. Our predictions show qualitative agreement with the data from [Hodby et al. , Phys. 
Rev. Lett. 94, f20402 (2005)] without the use of adjustable parameters. 



By ramping a magnetic field across a Feshbach reso- 
nance, loosely bound diatomic molecules have been cre- 
ated in two-component Fermi gases above 0, |2| and be- 
low ^ |E| the superfluid transition temperature. In 

] this process, there appear to be two distinct relaxation 
timescales: a "two-body" timescale in the range of lO's 
to lOO's of microseconds during which molecules form, 
and a longer "many-body" timescale in the range of lO's 
of milliseconds during which a condensate forms 0, |^ . 
A critical unresolved riddle emerging from these experi- 
ments is that the short-time molecular formation physics 
can not be fully accounted for by theories that treat the 
Feshbach ramp as a purely two-body process M, SK sug - 

[ gesting a many-body treatment is warranted |3, llO. HH ■ 
In this letter, we argue that in the above-cited exper- 
iments, multiple collisions per atom occur during the 
ramp. Thus, Feshbach molecule formation in a thermal 
gas must be viewed as a coUisional relaxation process, 

[ which we describe using kinetic theory 0, and ther- 

■ modynamics [Hill El. 

Several key experimental quantities in need of a com- 
prehensive theoretical understanding are illustrated qual- 
itatively in Fig. 1. As the magnetic field is tuned across 
the resonance, the molecular fraction x = ZA^M/A'tot in- 
creases to a maximum value x b ' '^h^re A^tot and Nm are 
the total and molecular populations. This behavior is 
represented in Fig. la by the error function Xerf(^) = 
Xb{1 " erf[-\/2(i3 — Bcen)l5Ey\/2, which has been used 
to fit experimental data 0, Fig. lb shows how the 
production of molecules increases and then saturates to 
a value less than unity as the ramp rate is decreased. 
An exponential function = 6xp(— Q;/|i?|)] gives 

a reasonable fit to data U, |3, M , however, a power-law 
, form can not be ruled out [l^- The rate constant a 
determines whether the sweep is adiabatic and xo is the 
maximum production efficiency obtained in the adiabatic 
limit a/|-B| ^ 1. In this letter, we identify the root cause 
of the saturation effect and determine the scaling behav- 
ior of the measured properties Been, 5B, a and xo- 

A useful idealization that is commonly invoked to un- 
derstand Feshbach molecule formation is to scale the sys- 
tem size down to two isolated trapped atoms initially in 
their ground state. This description is amenable to a 
straight-forward application of the Landau-Zener (LZ) 
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FIG. 1: Illustrations of the molecular fraction versus (a) mag- 
netic field B and (b) sweep rate \B\. 



level crossing theory 7], for which the effect of the mag- 
netic field ramp is easy to conceptualize: the state of the 
atom pair tracks the lowest energy state during the ramp, 
eventually crossing threshold to become a true two-body 
bound state. In the adiabatic limit, the LZ model pre- 
dicts that at the end of the sweep the probability of find- 
ing the atom pair in a bound molecular state is unity, 
corresponding to a production efficiency xo = 1- For fi- 
nite temperature gases, this two-atom LZ model is not 
directly applicable, and thus one does not necessarily ex- 
pect that a unit conversion efficiency can be achieved. 

Nevertheless, is there a way to incorporate the two- 
atom LZ theory into a model that takes into account the 
presence of many atoms in a trapped gas? The key to 
answering this is to realize that in a dilute gas there is an 
upper bound on the ramp time iramp set by the average 
time between collisions Tcoi: if iramp 3> Tcoi, an individual 
atom may undergo multiple collisions during the ramp, 
and this physics is not taken into account by the LZ ap- 
proach. In the opposite limit tramp ^ t-coI, each atom 
will likely encounter at most a single collision partner - 
its nearest neighbor, so that the LZ theory can be applied 
to isolated pairs of atoms distributed throughout the gas. 
Models working in this type of "static gas" approxima- 
tion (SGA) predict a maximum production efficiency of 
Xo = 0.5 lOj. Although this efficiency limit is con- 
sistent with the first two Feshbach sweep experiments in 
a Fermi gas 0, the more extensive recent study of 
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Hodby et al. measured values of xo in the range 

0.1 — 0.9, in disagreement with the SGA prediction. 

We can understand this apparent breakdown of the 
SGA by estimating Tcoi for typical experimental param- 
eters using T^~J — navy-ci, where a is the collisional 
cross section and Wroi is the average relative velocity. 
In the unitarity limit, a = Ana? /(I + k'^a?), which 
can be approximated by a = An/kp close to the reso- 
nance where the scattering length a diverges. Here kp 
is the Fermi wave vector. Estimating the density as 
n « {2mkBTp Z??)^^^ /6tt^ , where m is the atomic mass, 
and taking Vrei — A^Jk^TY /Tim^ the collisional time be- 
comes Tcoi ~ h/k^TY jS]- In the JILA experiments [l^ . 
typical Fermi temperatures are T-p = 350 nK, giving 
Tcoi = 20 ^s. With iramp 50 /iS in these experiments, 
multiple collisions per atom occur during the ramp and 
an alternative to the SGA is necessary. 

To go beyond the SGA, Hodby et al. devised a stochas- 
tic phase space sampling (SPSS) model based on the as- 
sumption that the probability of two atoms forming a 
molecule depends on their proximity in phase space 'IT"] . 
With a single fitting parameter - a cut-off radius in phase 
space - the phenomenological SPSS model predicts that 
that Xo is a universal function of poi that agrees remark- 
ably well with their experimental data. In this letter, 
we provide a theoretical backdrop for understanding this 
behavior and derive from first principles the universal 
functional form for Xo(poi)- 

Rather than try to adapt the two-atom LZ theory to 
the case of a finite-temperature gas, we have developed 



an entirely different many-body approach that is appro- 
priate for the regime tramp '3> Tco\- Using the Keldysh 
nonequilibrium Green's function formalism |l8j| . we de- 
rived a pair of coupled Boltzmann equations describing 
the dynamical evolution of the atomic f a (Pi t) and 
molecular fMiv^'i^Tt) distribution functions [l3|; we as- 
sume an equal spin mixture fA = f]=fi- Working 
within the quasi-particle approximation, we then devel- 
oped idealized models to study the dynamics and 
thermodynamics [Id^ of molecule formation. In this let- 
ter, we describe the salient features of this approach and 
apply our theory to the experiments of Hodby et al. |0| . 

In Tl^, we used our kinetic theory to derive rate equa- 
tions for the atomic (A) and molecular (M) populations 
Nj {t) = Jdr f dpfj (p, r, t)/h^, where j = A or M: 

NM{t) = lf{t)NA{t) ~ ld{t)NM{t), (1) 

and NA{t) = -Nuit). The quantity 7d(t) 
is the molecular dissociation rate and 7/(t) = 
J drnA(j^,t)jf{r,t)/NA{t) is the density-weighted aver- 
age rate of molecule formation. The density is defined by 
nj{r,t) — J dpfj{p,r,t)/h^. The rates of molecule for- 
mation 7/(r, t) and dissociation jdit) follow directly from 
the "in" and "out" terms of the atom-molecule collision 
integrals in the Boltzmann equations, which describe the 
relaxation toward chemical equilibrium due to resonant 
two-body collisions. Within the classical gas approxima- 
tion (i.e. fi <C 1), they are 
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y"dpi J dp2 y"dp3^(Pi +P2 - P3)'5(eA(Pi) + eA(P2) - ej\/(P3))/A(Pi)/A(P2), (3) 



where the explicit position and time dependence has been 
suppressed in the integrals. Neglecting self-energy ef- 
fects, the atom and molecule energies are e^(p, r) = 
p^/2m + [/cxt(r) and eA/(p,r,t) = p^/Am + 2J7oxt(r) + 
^rcs[B{t)], respectively, where [/oxt(r) is the external trap- 
ping potential and ercs[B{t)] is the renormalized energy 
of the resonant closed-channel state. This energy cor- 
responds to the peak position of the molecule spec- 
tral function with the above-threshold limiting forms 
ercs(S) = Afi{B — Bq) far from resonance \B — Bo\ ^ 
|A_B| and e-[cs{B) « t?/ma'^{B) close to the resonance 
\B — Bo\ ^ |Ai?| 25]. Here Bq is the resonance posi- 
tion, AB is the resonance width, and A/i is the mag- 
netic moment difference between the open and closed 
channels. The scattering length near a Feshbach reso- 
nance is a{B) ~ abg[l — AB/{B — Bo)], where Cbg is 
the background scattering length. Neglecting the energy 



dependence of interactions here, the Feshbach coupling 
strength g is related to physical parameters according to 
^2 = A'Kh'^ayjgAfj.AB/m. 



The rates in 
itly when the 
exp{-[ej(p,r) - 



and can be calculated explic- 
system is in equilibrium /j(p, r) = 
- iij]/kBT}, with ^iM = 2^iA = 2/i. 



The molecular dissociation rate then simplifies to the 
position-independent form 



Id 



2to1/2 

2 — ahgAiiABy/erc 



(4) 



Here, the unit step function ^(eros) follows from energy 
and momentum conservation and signifies that below 
threshold eres < 0, the pairs become stable and no longer 
can decay. We note that Q is consistent with the decay 
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rate derived in Ref. '19]. The formation rate reduces to 
7/(r) = 23/2p(r)e-^--/'=BT^^^ 

where p(r) = A^j^n^ (r) is the phase space density of the 
atomic component and Ath = /i/ V^TrmfceT is the thermal 
deBroghe wavelength. The density-weighted average is 
7^ = 7j^(0)/2'^/^. The formation rate vanishes as ei-cs 
is tuned below threshold, in the same way that 7^ does. 
This causes the saturation of molecule production during 
a Feshbach sweep: once e^es is tuned below threshold, no 
more molecules are formed via two-body collisions. 

Despite its simple intuitive form, Eq. cannot be 
solved without knowing the fj{p,r,t) of Eq. (O and 
Eq. Nonetheless, we can extract from Eq. an 

approximate form for the characteristic rate constant 
a that demarcates the adiabatic regime. At the very 
onset of molecule formation when x !> the solu- 
tion of Eq. ^ can be approximated by A''M(iramp) ~ 
(iVtot/2){l - exp[-/g''""''7/(r)dr]}. Making a change 
of variables from r to Cres in the integral, and taking 
the large detuning form e^cs = "I^mII-^I? obtain 
Nm{B) oc [1 — exp(— where the rate constant is 
a = /(J (iei.cs7/(ercs)/|A/i|. Evaluating the integral using 
the equilibrium form of 7/ gives 

a = 2^/^7r^^^ABn^o, (6) 

TO I A/i| 

where uao = 't-aIi" = 0) is the peak atomic density. This 
result l|njl resembles the Landau-Zener rate constant de- 
rived for a zero ternperature Bose gas in [^J. The ex- 
perimental data in |17| seem to be consistent with this 
linear dependence on density. 

In the adiabatic regime Q!/|i3| ^ 1, the system fol- 
lows an isentropic path as Crcs is ramped downward [lll| . 
When the resonance energy €rcs[B{t)] is ramped to neg- 
ative values, the two-body rates of molecule formation 
7/ and dissociation 7^ vanish, at which point the pro- 
cess of molecule formation must proceed via three-body 
collisional relaxation, which may play a role for very 
slow ramps j23|. When the timescale for three-body 
relaxation is much longer than iramp, it can be ne- 
glected Thus, molecule formation effectively 
halts when Cres crosses threshold at B = Bq, giving rise 
to the observed saturation of molecule production. 

In we calculated constant entropy contours for 
a quantum ideal gas mixture of fermionic atoms and 
bosonic molecules without accounting for this saturation 
effect. It is straightforward to incorporate this effect into 
our earlier thermodynamic calculation. For positive de- 
tunings where Cios > 0, the results are unchanged: the 
molecular fraction xi^rcs) is determined according to the 
conservation of entropy and total atom number with the 
constraint that the atomic and molecular components are 
in thermal Tj^i — Ta = T and chemical fiM = '^fJ-A = 
equilibrium. For negative detunings, the molecular frac- 
tion is held fixed x(ei.cs < 0) = x(0) = XO: where xo is 



0.4 




magnetic field {B - Bq) (10""* T) 

FIG. 2: Molecular formation curve for the case of '"'K with 
Ti/Tp = 0.5. The solid blue line in both graphs corresponds 
to the constant entropy equilibrium solution. The dashed 
red line Xerf (-B) with SB = 0.19 x 10"" T and Bcon ~ Bo = 
0.34 X 10"* T. The dashed blue line in the inset shows the case 
where the system remains in chemical equilibrium for below 
threshold. 



the maximum production efficiency. The details of this 
type of calculation are given in 'l5l|. 

In Fig. 2 we plot the molecular fraction along an isen- 
tropic path (solid blue lines) for the specific case of an 
initial temperature T-jT-p — 0.5, taking values for ^"K: 
So = 202 X 10~4 T, AS = 7.4 x IQ-* T, and Obg 
174oo 0. The inset shows x(£res) versus Crcs compared 
to the case where the gas remains in chemical equilibrium 
for negative detunings (dashed blue line). To plot x ver- 
sus magnetic field, the small detuning limit was used for 
ercs(^)- The resulting molecular formation curve x(^) is 
fit remarkably well by the error function Xerf (S) (dashed 
red line). The fitted width and centroid for four differ- 
ent initial temperatures T/T-p = {0.25,0.5,0.75, 1.0} are 
6B = {0.14,0.19,0.21,0.24} x lO""* T and Bcc„ - = 
{0.37, 0.34, 0.35, 0.36} x 10"'' T. These widths 5B are con- 
sistent with the value 0.2 x lO"** T measured in experi- 
ments 0,01, which is much smaller than the Feshbach res- 
onance itself (AS = 7.4 x 10"''' T). A crucial point arising 
from our analysis is that the centroid Been of the error 
function does not coincide with the resonance position 
Bo- The centroid occurs roughly where e-ccs{B) ~ kBTp, 
so that Been — Bo « ABahg^/mkBTF /h. 

Our saturation mechanism is quite general and occurs 
in both Bose and Fermi gases; examples of both sys- 
tems have been studied by Hodby et al. 17]. Within 
the classical gas approximation it is straightforward to 
relate xo to the initial peak atomic phase space den- 
sity poi. The maximum molecular fraction is given 
by Xo = 2NMiTf,pQf,ercs = 0)/A^tot, where T/ and 
PQf are the temperature and peak atomic phase space 
density at e^cs = 0. The total atomic population is 
.^tot = V^A + "^Nm, where for generality we allow for 
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FIG. 3: Comparison between theory and experiment. The 
solid blue curve is for a two-component Fermi gas including 
quantum statistical effects. The dashed green line shows the 
classical result, Eq. (|7J, for rj = 2. The solid black circles are 
experimental data for ''"K from Hodby et al. [l^. The two 
earliest experiments are those of JILA (red square) Q and 
Rice (red diamond) 0. The inset compares the ri = 1 case 
of Eq. (|7J (dashed green line) to data for **^Rb taken in the 
classical regime of low phase space density . 



the number of atomic spin components to be rj e {1,2}. 
The atomic and molecular populations can be written 
as Na{T,po) = {kBT/Tiufpo and Nm{T, po,e,^s) = 
{k-gT /hujY p'^eyiYi{—e^as/kBT), where the mean trapping 
frequency is (D = {uJxi^y'jJzY^'^ ■ Using these expres- 
sions, the molecular fraction takes the form xo = 
Pofli'nl'^ + Po/)- We can relate po/ to poi using en- 
tropy conservation. The atomic and molecular entropies 
can be written as S'A(r, po) = fcB-^yi(4 — Inpo) and 



SM{T,po,eres) = kuNMii + eres/ksT - 2\npQ). Sub- 
stituting these into the entropy conservation equation 
riSA{Ti,poi) = r]SA{Tf,pof)+SMiTf,pof,e,.cs = 0) yields 
the relation In pAoi = In pAof + 2xo- From these observa- 
tions, the following key result can be derived 

2Xo + lnf^-^) =lnpo.. (7) 
V2 1 - xoj 

The solution xo of this transcendental equation is a 
universal function of the the initial peak atomic phase 
space density poi aud the number of spin components rj. 
A similar derivation for an ideal quantum gas results in 
a much more complicated expression, in which xo a-lso 
depends on the quantum statistics of the different com- 
ponents. In Fig. 3 we compare our theoretical prediction 
of Xo, both with (solid blue) and without (dashed green) 
quantum statistical effects, to experimental data and find 
qualitative agreement with the JILA dataQ, 0| but not 
with the Rice data point (red diamond) . 

In summary, we have treated Feshbach molecule for- 
mation in a dilute gas during a magnetic field ramp as 
a relaxation process described by kinetic theory [l^, [131 
and thermodynamics . Such a treatment is appropri- 
ate to the conditions of present experiments 1, 2, l?!. 
Our key findings are: (i) the rate factor a is proportional 
to density, (ii) the width 6B and centroid Been of the 
formation curve are related to the temperature T of the 
gas, and (iii) the molecule production efficiency xo is a 
universal function of the initial peak atomic phase space 
density poi i which in the classical regime is approximated 
by the ideal gas result Eq. {Tj). Our predictions are con- 
sistent with most of the existing experimental data. 

We thank E. Hodby and C. Regal for providing the 
experimental data used in Fig. 3. We appreciate useful 
discussions with T. Nikuni, P. Julienne and E. Tiesinga. 
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